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We study maximal bipartite entanglement in valence-bond states and show that the average bi- 
partite entanglement E^, between a sub-system of two spins and the rest of the system, can be 
maximized through a homogenized superposition of the valence-bond states. Our derived maxi- 
mal ijjj rapidly increases with system size and saturates at its maximum allowed value. We also 
demonstrate that our maximal states are ground states of an infinite range Heisenberg model 
(IRHM) and represent a new class of resonating-valence-bond (RVB) states. The entangled RVB 
states produced from our IRHM are robust against interaction of spins with both local and global 
phonons and represent a new class of decoherence free states. 

PACS numbers: 03.67.Mn, 03.65.Yz, 03.67.Bg, 05.50.-|-q 
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I. INTRODUCTION. 

Quantum entanglement, a manifestation of non- 
locality, is a precious resource for quantum computation 
and quantum information [l[ and signifies correlations 
in many-body systems. Quantum algorithms that would 
significantly accelerate a classical computation must rely 
on highly entangled states since slightly entangled states 
can be simulated efficiently on a classical computer 0- 
The strength of correlations of fluctuations of observ- 
ables (such as density, magnetization, etc.) in a many- 
body system is a reflection of the degree of entanglement 
(for pure states) 0]. Thus, characterization of multi- 
particle entanglement and production of maximal/high 
multi-qubit entanglement is vital for quantum compu- 
tational studies and for mutual enrichment of quantum 
information and many-body condensed matter physics. 

Intensive work on entanglement during the past decade 
has led to the proposal of numerous measures of entan- 
glement [3, [Hi- While two-party entanglement is quite 
well understood, entanglement in a multi-party system 
is an area of immense current interest |6l-ll3{. In the 
quest for maximally entangled states, so far only few- 
qubit maximally entangled states such as two-qubit Bell 
states, thrce-qubit Greenberger-Horne-Zeilinger (GHZ) 
states, and four-qubit Higuchi-Sudbery (HS) states [ij 
have been clearly identified. It would be of considerable 
interest to generate these maximally entangled states as 
the ground states of a physically realizable spin model. 
While the GHZ states could be obtained as the ground 
state of an anisotropic Heisenberg model [l5| , maximally 
entangled four-qubit and five-qubit states could be ob- 
tained only as excited states jla ]. 

Decoherence is one of the main obstacles for the 
preparation, observation, and implementation of 
multi-qubit entangled states. Since coupling to the 
environment and the concomitant entanglement fragility 
are ubiquitous [l], it is imperative that progress 

be made in understanding decoherence as well. In 



the past decoherence-free-subspace (DFS) [H, [T^ has 
been shown to exist in the Hilbert space of a model 
where all qubits of the quantum system are coupled 
to a common environment with equal strength. This 
situation manifests when the distance between the 
qubits is negligible compared to the correlation length 
of the environment. If the dynamic symmetry of the 
system-environment interaction selects a set of orthonor- 
mal basis vectors (of the reduced Hilbert space) that 
is unaffected by environmental interaction, then such a 
subspace is called a DFS. An important application of 
decoherence free subspaces lies in developing quantum 
error correcting codes [13, HH . These subspaces prevent 
the loss of information due to destructive environmental 
interactions and thus circumvent the need for stabiliza- 
tion methods for quantum computation and quantum 
information. Decoherence free states have also been 
shown to be useful for quantum communication be- 
tween parties without a common reference frame (22| - [26t . 



RVB states have provided interesting insights for 
understanding strongly correlate d p henomena such as 
physics of high Tj;, cuprates [13, [11], superconductivity 
in organic solids [29[ , insulator-superconductor transition 
in boron-doped diamond [lOl, etc. Furthermore, RVB 
states have also been proposed as robust basis states 
for topological quantum computation [3l|. The multi- 
particle entanglement has also been investigated in RVB 
states that were proposed as states close in energ y t o 
the ground state of the Heisenberg Hamiltonian |32h34| . 
However, there has been no explicit construction of RVB 
states that would represent maximally entangled valence 
bond states. In this paper we construct a new class of 
RVB states that are ground states of IRHM, that have 
high E'^ entanglement, and are decoherence free. Our 
IRHM Hamiltonian couples to an environment that dis- 
tinguishes between the qubits unlike the case considered 
in Ref. [l^ We believe that an improved understanding 
of the entanglement and decoherence properties of RVB 
states will enable their implementation for quantum com- 
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putation and quantum information purposes. 

The remainder of this paper is organized as follows. 
In section II, we introduce entanglement entropy and 
demonstrate explicitly that isotropy and homogeneity 
maximizes entanglement between two spins and the 
rest of the system in the two limits of even-N-qubits, 
i.e., for N = 4 and N ^ oo. For intermediate cvcn- 
N-spin states, our proposed entangled states are only 
shown to maximize E^, among valence-bond states. In 
section III, from the ground states of our IRHM, we con- 
struct explicitly states with maximal E'^ entanglement 
among valence-bond states and show that these states 
form a special class of RVB states. Next, in section IV, 
we analyze the robustness of this special class of entan- 
gled states. We show that our IRHM Hamiltonian, even 
upon inclusion of both local and global optical phonons 
that are coupled to spins, produces ground states that 
are decoherence free. Lastly, in the final section V, we 
conclude after commenting on the physical realizability 
of our proposed high E^, entangled states. 



II. ENTANGLEMENT FOR TWO-QUBIT 
REDUCED DENSITY MATRIX IN ISOTROPIC 
SYSTEM. 

For a bipartite system AB in a pure state, von Neu- 
mann entropy Ey measures the entanglement between 
the subsystems A and B. From the reduced density ma- 
trices Pa = trsp^^ and ps = trAP^^ obtained from the 
pure state p'^^, we obtain 

Ev{pa) = -tr{pAlog2 pa) 

= -tr{pB log2 pb) = E.„{pb). (1) 

Using the basis | |), | t), and = icr', the single qubit 
density matrix can be written as [35j 
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(2) 



Throughout this paper, we consider only states |^') that 
are eigenstates of the z-component of the total spin oper- 
ator {S^^^g^l) with eigenvalue 5|.; furthermore, we focus 
on only isotropic states. It then follows that, {S^) = 
and {Sf) = leading to the single-qubit density matrix 
to be maximally mixed and thus maximizing entangle- 
ment Ey{pi). On realizing that (Sf) = 0, we obtain the 
following expression for the two-qubit reduced density 
matrix 1351: 
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(3) 

Here, isotropy implies 0.5(S',"5+) = 0.5{S+S') = 



{S^S^). Thus, the von Neumann 



entropy Ey [pij ) can be expressed as 



1 



Ey{p,,) = 2 - - [3(1 + 4(5f log2(l + 4(5f 

+ (1 - 12(Sf log2(l - l2{StS^^))m 
which implies that 



For our states, S^^i^i\^) 
(5fE,^|> =0,thatis. 



(5j Sj 



-{St 



(5) 



We will now maximize the total entanglement entropy 
Till j^i Ev{pij) subject to the above constraint in Eq. ([5|). 
To this end, we will employ the method of Lagrange mul- 
tipliers and define the Lagrange function A as follows: 



K=Y,Ey{p,,)-Y,\AY,{sts^^ 



Then, setting 
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(6) 



dA 



yields 



A; = 3 loe 



(1 



UiSfSt, 



)) 



(l + 4(5f5^)) 



(7) 



which implies that the optimal (SfSj) is independent of 
J for j ^ i. Consequently, it follows from Eq. ([5]) that 
J2i,j^iEv{pij) is maximized when = -j^^ziYj, 

i.e., when the isotropic state is homogeneous. The av- 
erage entanglement entropy between the subsystem of 
two spins and the rest of the system (of — 2 spins) 
E^ = [1/A''(A''- 1)] J2i.j^i EviPij) has a maximum value 
given by 



E. 
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^°S^li-4(A^ 
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log2 T 



4(A - 1) 



(8) 



It is interesting to note that for A oo, the above 
expression yields E^ ^ 2. In fact, E'^ approaches the 
maximum possible value 2 quite rapidly as can be seen 
from Fig. [TJ Next, we observe that for A = 4, our ex- 
pression for i^^^TOa^. in Eq. ^ yields the same entangle- 
ment entropy value 1-1-0.5 log2 3 as that obtained for the 
four-qubit maximally entangled HS state [l^, ■ Fur- 
thermore, our approach explains why all the pairs in the 
HS state give the same entanglement value. Contrast- 
ingly, for the isotropic ground state, in the case of a 
Heisenberg chain with nearest-neighbor interaction, for 
four spins the entanglement entropy E'^ ~ 1.21 (with 
(SfSj) = —0.5/3) while for an infinite chain the von Neu- 
mann entropy E^ = 1.37 (with (SfS^) w -0.443/3) both 
of which are far less than our maximal E'^ values above. 
It is of interest to note that, when A = 4 or A — ^ oo, 
the maximal values of E'^ for isotropic systems are the 
same as the maximal E'^ values for general systems (i.e., 
systems that can be either isotropic or non-isotropic). 
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FIG. 1. Normalized entanglement E^/Eniaxj for two-qubit 
reduced density matrix, measured by (a) von Neumann en- 
tropy (Ev/Ev,niax) and (b) i-concurrence (Ic/Ic,max), for N- 
qubit valence-bond systems. 



We also note that homogeneity, under the constraint 
of Eq. ([5]), maximizes the i-concurrence (/c) [36j given 
by 



NiN - 1) ^ 

i,j>i 



^2[l -Tr{pl)]. (9) 



As shown in Fig. 1, Ic also monotonically increases with 
system size. Although we considered von Neumann en- 
tropy and i-concurrence as entanglement measures, our 
homogeneous states should also maximize other measures 
of entanglement for pij in valence-bond systems. 



entanglement is given by Eq. It is important to note 
that, while all possible St = states are ground states of 
IRHM, for some Hamiltonians only some of the possible 
St = states are ground states. For instance, for the 
nearest-neighbor antiferromagnetic Heisenberg Hamilto- 
nian or for the Majumdar- Ghosh model Hamiltonian (STj 
only some of the possible St = states are ground states 
and consequently (for these systems) it is not possible to 
construct a homogenized linear combination of St = 
ground states that produces maximal possible entangle- 
ment El„,^^ given by Eq. ©. 

We will now proceed to construct entangled states for 
spins that maximize E^. We first note that there are 
(A — 1)!! eigcnstates with 5*^ = for the Hamiltonian of 
Eq. ([ini). Each of the St = states is a product of N/2 
two-spin singlet states of the form | fi) — | it) (with no 
pair of singlets sharing a spin). Of these (A — 1)!! prod- 
uct combinations with St = 0, only ^Cn_ — ^ C n__i = 

A!/[(A/2)!(A/2 + l)!] products are linearly independent. 
A particular set of linearly independent St = states 
are the Rumer states [s^ [s^ which are made up of non- 
crossing singlets. For A^ = 6, the Rumer states are the 5 
diagrams with non-crossing singlets shown in Fig. 2. 
For TV = 4 spins, two linearly independent St = 
f2'"°) ® l*ff and |$ff=°) ^ 

S ■ ■ —0 

<I>jj-^ ) is a two-spin singlet state for 
spins at sites i and j with Sij being the total spin of Si 
and Sj. It is worth noting that 



eigenstates are 1$?,^""°) ® 1$?^"°) and |$fi 
1^23'"°) where 



l'i>ff=°> 



(12) 



III. MULTI-QUBIT ENTANGLED STATES. 



A. Highly entangled ground states of IRHM 



Using the above relation one can establish that there are 

only ^Cjv — ^Cjv linearly independent St = states. 

2 2 

Furthermore, we also observe that 



[S1.S3 + S2.Si + Si.Si + Ss.Sall^fs' 



Sl2=0\ 
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S34=0\ 
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We will now demonstrate that ground states of the 
IRHM Hamiltonian 



J 



IRHM 



=jj:s^s,^^ 



i,j>i 



,(10) 



will produce the same amount of entanglement as given 
by Eq. ([8]). We note that -ffiRHM commutes with both 

S^Totai and (E. ^1)' ^lotai)- Eq. m, it is under- 
stood that J ~ J* / [N ~\) (where J* is a finite quantity) 
so that the energy per site remains finite as A^ — > 00. The 
eigenstates of i?iRHM correspond to eigenenergies 



E. 



St 



St {Sj 



1) 



3A^ 



(11) 



where St is the total spin eigenvalue. The ground state 
corresponds to St = which is rotationally invariant 
and also implies that S^ = 0. Thus for a homogenized 
linear combination of the St = states of i^iRHM, the 



using which it can be shown that any St = state made 
up of A/2 singlets will be an eigenstate of Eq. pO)) . 

Now, while the two symmetry broken St = singlet 
states [on the right hand side of Eq. ([T2l) ] are isotropic, 
they are not homogeneous. To make our entangled state 
homogeneous, we construct the following state that has 
the symmetry that {S^' Sj) is the same (= —1/12) for all 
pairs i and j and any rj (= x, y, z): 



l<i> 



S34=0 

34 



»+^3(l* 



Sl4=0x 

14 / 



23 



^[(inn) + unt)) + ^3(imt) 



;m)) 



-^3 (im;) + u;n))] 



I* 



Hs; 



(13) 



where ws's are cube roots of unity. The above |4'hs) 
state is the maximally entangled four-qubit HS state 
studied by many authors [lol Tisl [T^ . [l6| . However, it 
has not been shown earlier that I'I'hs) can be expressed 
as a superposition of St = states. Importantly, the 
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FIG. 2. Homogenized linear combination of the St = states 
in both (a) and (b) give the same maximal entanglement for 
six qubits. 



model. 

The above strategy of constructing homogeneous and 
maximal states from isotropic four-qubit states can be 
extended to the six-qubit case as well using the St — 
ground states of -//irhm shown in the six diagrams of 
Fig. [21 The entanglement E'^ for the six-qubit states is 
maximized by taking suitable resonance hybrids of the 
diagrams shown in Figs. [5] (a) or (b). These resonance 
states homogenize the ground state in the sense men- 
tioned above (i.e., produce site independent correlation 
hmctions (S'-'5'') = -1/20): 
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(14) 



and 



half-filled large U/t Hubbard model on a regular tetra- 
hedron yields the above two broken-symmetry St = 
states as ground states from which the HS-state can be 
constructed. In previous reported works [i^ |4l|, al- 
though the authors studied entanglement in the model 
H = Jj2i j=i 2 3 4 Si-Sj, they did not point out that the 
maximally entangled HS state is a ground state of the 



= W4( 


I*f2" 




5 |$3t 




5 I'J'fr 


=°)) 
















+ w|( 






5 |ci>ff 




5 |<i>ff 





1*^ 



which can be rewritten as 
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mm) - 1 mm)) + ^4{\ mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) - 1 mm)) 



-Loiii mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) - 1 mm)) 
■ w|(i mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) - 1 mm))], 
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mm) - 1 mm)) + ^4(1 mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) - 1 mm)) 



■ c^Kl mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) - 1 mm)) 

■ ^1(1 mm) + 1 mm) + 1 mm) - 1 mm) - 1 mm) - 1 mm))] 



(17) 



where 0.4 's are fourth roots of unity. The von Neumann 
entropy E^,, for two-qubit reduced density matrix ob- 
tained from l^-g), (l^-g))*, l^-^), or (j^-^))*, is 1.921964 
which is the same as iJ^^maa; proposed by our general for- 
mula in Eq. ([5]). The conjecture made by Brown et al. 
[l3| that the multi-qubit maximally entangled states al- 
ways have their reduced single qubit density matrix max- 
imally mixed is satisfied by our states J^I'e)- It should be 
noted that the highly entangled six-qubit state reported 
in Ref. [lo| . although yields a higher entanglement, is 
not an eigenstate of the S^^f^i operator. 

Cabello [l^ has constructed supersinglcts of four and 
six qubits which are decoherence free. It is interesting 
to note that they form the ground states of our IRHM 
Hamiltonian. However, the average two particle von Neu- 



mann entanglement entropy E'^ for these supersinglets 
is less than the entanglement for our states constructed 
above. In the case of the four qubit supersinglet state 

iS^^''^, the entanglement entropies for all possible bi- 

partitions arc given by Ey(pi2) = Ey^p^^) = 1.5849, 
Evipis) = Ey{p24) = 1.20 75, Ey{pu) = ^.(^23) = 
1.2075 which together yield the average two particle en- 
tropy value 1.3333, i.e., a quantity clearly smaller than 
the von Neumann entropy E'^ = 1 + 0.5 loga 3 sa 1.79248 
obtained for the four-qubit maximally entangled HS 

the 



state. 



c(2)\ 

Dg ) case, 



In the six qubit supersinglet 

various bipartitions produce the average entropy E^ = 
1.657996 which is noticeably smaller than the entropy 



value E:, 



~ 1.921964 obtained for our six qubit en- 
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tangled states. 

On the issue of how many maximal E'^ states are pos- 
sible for N qubit systems, we would like to say that it 
cannot exceed the total number of linearly independent 
St = states (i.e., the total number of Rumer states ). 
As regards existence of homogeneous states for > 8 
qubit systems, wc provide arguments in Appendix A. 

B. Resonating- valence-bond picture 

The construction of our entangled maximal E'^ states 
from the ground states of IRHM can be visualized by 
using a RVB picture. Our maximal E^ states can be 
regarded as a new class of RVB states made of homog- 
enized superposition of isotropic St = valence bond 
states. We will now compare the entanglement proper- 
ties of our RVB states and the general RVB states |^')rvb 
of Ref. m given below: 

|^')rvb ^ f{il,-;iM,jl,-;jM)\{il,jl)--{iM,jM)), 

where M represents the number of sites in each sub- 
lattice and / is assumed to be isotropic over the lattice. 
Also, liikjk)) = 75(l5)iJ -^k- I- denotes 
the singlet dimer connecting a site in sub-lattice A with 
a site in sub-lattice B. The valence bond basis (used for 
the above RVB state |^')rvb) form an over complete set 
while our RVB states are constructed from a complete 
set of^Cw — ^Civ_i states. 

The rotational invariance of the two qubit reduced den- 
sity matrix of the RVB states allows us to write them in 
the form of a Werner state: 

p^{p)-p\mm\ + ^h, (18) 

where for 1/3 < p < 1 the Werner state has the spins 
at i and j entangled with each other. For the special 
case of the "RVB gas" for the |4')rvb state, where / is 
a constant (corresponding to equal amplitude superpo- 
sition of all bipartite valence bond coverings), one gets 
the exact result p = 1/3-1- 2/(3A^) which implies that 
all finite size systems have a non-zero tangle (or en- 
tanglement) between the two sites [1^. Next, for the 
"RVB liquid" case (involving equal amplitude superpo- 
sition of all nearest-neighbor singlet valence bond cover- 
ings of a lattice), Monte Carlo calculations extrapolated 
to the thermodynamic limit yield p = 0.3946(3) > 1/ 3 
[33| . i.e., a non-zero tangle between the two sites [43| . 
In contrast, our maximal E'^ RVB states yield zero en- 
tanglement between the two spins for all system sizes as 
demonstrated below. It has been shown that the SU (2) 
symmetry of the RVB states ensures that the two-spin 
correlation function and the parameter p of the Werner 
state are related as 

{<i>\S..S,\^) = -^p. (19) 



Then, since our RVB states produce 

{s^s^) = {sfsy) = {s^s;) = ^^~[^y (20) 

it follows that p = and thus for all even A^ > 4 

systems wc get zero ent ang lement between the two sites. 
Lastly, based on Ref. [Sj, wc find that the monogamy 
argument yields the bound p < 1/3-1- 2/ {3\/N — 1) 
while the quantum telecloning argument produces the 
bound p < 1/3 + 2/[3(A^ — 1)]. Compared to our exact 
value of p = 1/(A^ ^ l)i both these bounds are weaker 
bounds (i.e., show zero two-site entanglement with cer- 
tainty only in the thermodynamic limit). Thus (among 
various RVB states) wc sec that our entangled high E^ 
RVB states, while producing maximum entanglement be- 
tween a pair and the rest of the system, yield zero entan- 
glement among the two spins of the pair. 

We will now remark on the entanglement of any set of n 
(> 2) sites with the rest of the system. The valence bond 
entanglement entropy picture of Alet et al. in Ref. [13 
uses the number of valence bond two-spin singlets shared 
by the two subsystems as a measure of their entangle- 
ment with each other. Thus, our 'homogenized' RVB 
states, will always have shared bonds between the two 
subsystems and thus show high bipartite entanglement. 
However, exact quantitative analysis needs to be carried 
out. Extending our approach to deriving the n-qubit re- 
duced density matrix in terms of n-particle correlation 
functions (when n > 2) and obtaining entanglement ex- 
pressions and wavcfunctions corresponding to maximal 
bipartite entanglement is a non-trivial exercise and is left 
for future studies. 



IV. ROBUST ENTANGLEMENT. 

The real quantum computer will not be free from noise 
and thus the entangled states have a tendency to undergo 
decoherence. However if the entangled state, in the pres- 
ence of potentially decoherence producing interactions, 
still remains entangled we call such a state robust oth- 
erwise it is fragile. We will now present two extreme 
cases of system-environment interaction scenarios where 
the ground states of our IRHM as well as the highly en- 
tangled RVB states (that we construct from the ground 
states) are decoherence free. 



A. Decoherence due to local optical phonons 

To study decoherence due to phonons, we consider 
interaction with optical phonons such as would be en- 
countered when considering a Hubbard model. The total 
Hamiltonian Ht is given by 

Ht = HiBSM + g^Q ^ S-{al + a,) + Wo ^ ala-i, (21) 

i i 
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where a is the phonon destruction operator, luq is the 
optical phonon frequency, and g is the couphng strength. 
Then, using the steps given in Appendix B, we obtain 
the following effective Hamiltonian He through second- 
order perturbation theory for strong coupling {g > 1) 
and non-adiabatic ( J/wq i£ 1) conditions: 



i 

(22) 



where 



2 Pe-'^y 
J^ = Je-'-> -{N-2)h{g)— , (23) 



J|| = J+[2A(5) + /2(5)] 



2wn 



(24) 



with h{g) = Er=i5'"/("!'^) and ^(.g) = 
Y.n=iYZ=i9^^^^"'V[n\m\{n + m)]. It is inter- 
esting to note that the eigenstates of the effective 
Hamiltonian He in Eq. ([22]) are identical to those of 
the original Hamiltonian i?iRHM in Eq. (jlOp because 



Eij>j('S'i''S'|),^^iRHM] = 0. Furthermore, even upon 



carrying out higher-order (i.e., beyond second-order) 
perturbation theory (as discussed in Appendix B) , we 
still get an effective Hamiltonian Hef f of the following 
form that has the same eigenstates as the IRHM. 



H. 



i,j>i 



(25) 



where Jxy and Jz are functions of the S^^^^i (= J2k ) 
operator. It is the infinite range of the Heisenberg model 
that enables the eigenstates of the system to remain un- 
changed. Furthermore, the ground state can also remain 
unchanged (with St — and = 0) if the effect of 
the term —C^^Sf with C > [such as the last term 
in Eq. ([H])] is canceled by a magnetic field. Thus the 
set of linearly independent ground states form the lowest 
energy subspace that is immune to decoherence. Further- 
more, our highly entangled RVB states constructed from 
these ground states [such as those given by Eq. ([T^ for 
four qubits and by Eqs. P^ - (fT7|) for six qubits] are also 
free of decoherence. Next, we study the decoherence in a 
dynamical context also and see how such ground states 
can remain robust and constitute a decoherence free sub- 
space. 



B. Dynamical evolution and DFS 

The robustness of entanglement in a system can also 
be looked from the dynamical perspective. We consider 



the following total Hamiltonian where all qubits of our 
IRHM interact identically with the environment. 



Hrot = ffiRHM+y^ ^ gkUJk{a\ + ak) 



E 



(26) 



The dynamics of the system can be studied through the 
following non-Markovian master equation for the reduced 
density operator ps{t) [i^ : 



dpsjt) 
dt 



.-t[Hs,ps{t)]+F{t)[Lps{t),L] 
+F*it)[L,psm], 



(27) 



where Hs is the Hamiltonian of the system and L is the 
system operator that couples with the bath and satisfies 
the constraint [L, Hg] = 0. For the total Hamiltonian in 
Eq. Hs = HmHM and L = J2^S^. Also, F{t) = 

/p a{t — s)ds where a{t — s) = T]{t — s) + iiy{t — s) is the 
bath correlation function at temperature T with 

V{t - s) = E ^^''^^ cos[ujk{t - s)], 

vit-s) = ~J2\gk\^sm[uk(t^s)]. (28) 



The function F{t) governs the non-Markovian dynamical 
features of the system. 

Let {|n)} be the eigen basis in which both the operators 
L and Hg are simultaneously diagonalized. Upon solving 
the master equation explicitly we get [4^ : 

Prnnit) ^{n\ps\m) 

=exp (-J[(i?„ - Ern)t + [ll - ll)Y{t)]) 

X exp[-(;„ - lrafX{t)]pran{()), (29) 

where En and In are defined through iJiRHMj^) — En\n) 
and X^j'^'i^l") — ln\n)- Furthermore, X{t) = FR{s)ds 
and Y{t) = /J Fi{s)ds with Ffl(t) + iFi{t) = F{t). This 
implies that, when states |m) and |n) have the same to- 
tal spin 5*7^ and the same z-component of the total spin 
( i.e.. En = Em and = Im), the matrix elements 
Pmn{t) = Pnm(O) and the states \m) and \n) belong to 
DFS. Furthermore, the density matrix for two qubits at 
sites i,j (obtained from the groundstates of our SU(2) 
symmetric IHRM Hamiltonian) exhibits time indepen- 
dence, that is, pij{t) = pij{0) with its elements given 
by Eq. This shows that the dynamics of the sys- 

tem has no effect on the entanglement of the ground 
state of the system which implies that these states are 
highly /maximally robust. Consequently, our high E^ 
RVB states constructed from the ground states of our 
IRHM are also free of decoherence. In future, using the 
above framework, we will consider the interesting case of 
dynamical evolution and decoherence of a superposition 
of states with different values. 



V. CONCLUSIONS. 

Although both spins in a two-spin singlet state | ti 
) — I 4,t) are monogamous (i.e., they cannot be entangled 
with any other spin), by using a homogenized superposi- 
tion of valenee bond states (eaeh of whieh is a product of 
N/2 two-spin singlets), we managed to distribute entan- 
glement efficiently such that any pair is maximally en- 
tangled with the rest of the RVB system while concomi- 
tantly the constituent spins of the pair are completely 
unentangled with each other. Thus we get the converse 
of monogamy for a pair of spins! Now, while total spin 
zero states are quite commonly ground states (as shown 
by Lieb-Mattis theorem [1^), it has not been recognized 
that one can generate high bipartite entanglement 
from such states. Our maximal E"^ RVB states are phys- 
ically realizable in systems such as infinite-range large 
U /t Hubbard model and infinite-range hard-core boson 
model with frustrated hopping [46| when they are at half- 
filling. 

Our high E'^ entanglement RVB states are free of de- 
coherence when the system interacts locally with opti- 
cal phonons. Furthermore, the set of orthogonal ground 
states of our IRHM form a DFS and evolve unitarily when 
all qubits are exposed to the same decoherening collec- 
tive noise. Thus our maximal E^ RVB states are free of 
decoherence for two extreme types of coupling with the 
environment! Our decoherence free RVB states can be 
used for error free quantum computation and quantum 
communication. 

In summary, our maximal E'^ RVB states have led us to 
IRHM for their realization and IRHM in turn generated 
maximal E^ states as highly robust ground states. 

Appendix A: Existence of homogeneous states for 
^(^8) qubit systems 

Let I '(/'„) be a set of linearly independent St = 
states with number of linearly independent states til — 
(Af/2)!(^/2-i-i)! • Then the homogeneous ground state can 
be written as: 

riL 

|*ifGs) -^(«n+i/3„)IV^n> (Al) 
n=l 

where a„ and /?„ are real. We have ^ C2 number of 
{SfS^) correlation functions appearing in the expression 
for El. For iV > 8, we note that < 271^; in fact, the 
ratio 2nL/^C2 is unity for TV = 8 while for iV > 8 the 
ratio monotonically and rapidly increases. In order to 
maximize the system E^., we must have {S^S^) indepen- 
dent of i and j and thus make the system homogeneous. 
Homogeneity condition implies that the equality of the 
^ C2 correlation functions (SfSj) produces ^6*2 — 1 in- 
dependent equations. Now, in \'^hgs), there are only 
271^ — 1 independent a„ and /3„ because of the normal- 
ization constraint. Furthermore for N > 8, in the ^C2 — 1 
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i j k 1 i j k 1 

(b) (c) 

i j k 1 i j k 1 

(d) (e) 

' '(f)' ' 

FIG. 3. Open loop hopping processes contributing to ef- 
fective hopping term T^* in third-order perturbation theory. 
Here empty circles correspond to sites with no particles while 
filled circles correspond to sites with hard-core-bosons. The 
numbers 1,2, and 3 indicate the order of hopping. 



equations involving (S'f S'|), the number of independent 
coefficients is greater than or equal the number of equa- 
tions, i.e., 2nL — 1 > ^C2 — 1. Thus there are enough 
number of unknowns to produce homogeneity. 

Appendix B: Third-order perturbation for IRHM 
system coupled to local phonons 

The total Hamiltonian Ht, involving interaction with 
local optical phonons, is given by 

Ht = i^iRHM + gi^o ^ Si{al + Qi) + Wo ^ aJai,(Bl) 

i i 

where a is the phonon destruction operator, ujq is the 
optical phonon frequency, and g is the coupling strength. 
Now. we make the connection that the spin operators can 
be expressed in terms of hard-core-boson (HCB) creation 
and destruction operators and b, i.e., b^ — 5*+, b = S~ , 
and b^b = -1-0.5. We then observe that conservation 
of 

^Totai implies conservation of total number of HCB. 
The total Hamiltonian is then given by 

H = [O.bbjbj + H.c. + n,nj] 

i,j>i 

+Wo^a]aj +5a;o^"j(aj +a]), (B2) 
j 3 

where Uj = b^jbj. Then we perform the well-known Lang- 
Firsov (LF) transformation [U 13 on this Hamilto- 
nian. Under the LF transformation given by He^^ = 



8 





(b) 



(b) 





2 
(c) 



FIG. 4. Closed-loop hopping processes contributing to efTec- 
tive interaction term in third-order perturbation theory. 
Here filled (empty) circles correspond to sites with (without) 
hard-core-bosons. The numbers 1, 2, and 3 represent hopping 
sequence. 



FIG. 5. Hopping processes (involving closed loops) contribut- 
ing to effective hopping term Tq„ in third-order perturbation 
theory. Filled (empty) circles represent occupied (unoccu- 
pied) sites. 



Hq + Hi with S = —g^i ni{ai — aj), the operators bj 
and Uj transform like fermions and phonons in the Hol- 
stein model. This is due to the interesting commutation 
properties of HCB given below: 



[b,,bj] = ==0, for I ^ J, 

{Kbl} = l. 



(B3) 



Next, the unperturbed Hamiltonian Hq is identified to 
be m 

Ho = [(0.5e-9'6|6j 4- H.c.) + n^n^] 



i.]>i 



(B4) 



For strong coupling [g > 1) and non-adiabatic (J/ujq < 



I) conditions, on noting that — Eq^ is a positive 
integral multiple of ujq and E!^"'' - E^"'' ~ Je-f' < ljq, 
we get the following second-order term [i^ 



where 



t{2) 



E 



[(0, 



-{N-2)h{g) 



2^0 



H.c.) + j\i^^n,nj 



-{N 



(B7) 



while the perturbation Hi is given by 

Hi = jY^ [^-^e-^'blbj + Y^.c.]{sfs'l - 1}, (B5) 



i,j>i 



where 



,)] and g'^ojQ is the polaronic 



exp[±g(ai 
binding energy. 

Since the unperturbed Hamiltonian Hq does not con- 
tain any interaction terms, we represent its eigenstates 
as \n,m) = |n)/ich ® \m)ph with the corresponding 
eigenenergies En,rn = E^n^ + EP^^ . On observiug that 
(0, 0|i7i|0, 0) = 0, we obtain the next relevant second- 
order perturbation term (48| 



(0,0|7Ji|n,TO)(n,TO|i/i|0,0) 



0,0 



E„ 



(B6) 



(2) 



[2/1(5) + /2 (.9)] 



2wn 



(B9) 



with /i(5) = Er=i5'"/("!«) and h{g) = 
Er=iEm=if^^""^'"V["!"i!(" + "^)]- The small pa- 
rameter for our perturbation theory is t/{guJo)- We now 
make the important observation that the eigenstates 
of the effective Hamiltonian Hq + H^-^^ are identical to 
those of the original Hamiltonian -ffiRHM in Eq. (|10|) 
because J2i,j>ii'^i'^j) and i/iRHM commute. 

Next, we will show that the third-order perturbation 
theory also produces a term that has the same eigen- 
states as IRHM. To this end, we obtain the following 
third-order perturbation term in the effective Hamilto- 
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The second hopping process Tj* m Fig. [3] (b) is given by 



ph{0\Hi\m)ph ph{m\Hi\n)ph ph{n\Hi\0)ph 



(BIO) 



Here AE„, = EP^ - Ef^^ . Evaluation of H^^^ leads to 
various hopping terms and interaction terms. 



6 3 



triTji T- / ^ i^cn-^ Cn 



n=l 



n=l 



(Bll) 



where t„ - (J^e 9^)/{g'^ujo)'^, ten ^ J^e s^/(gujo)^, and 
Vn ~ /(g'^^o)'^ (as will be explained later). We will 
demonstrate below that iJ*^^' is of the following form 



^(3) = j2 



T(^nfc)6|6,; + H.c. 



,l>i L k 



EnE-^-> 



(B12) 



where T and V are functions of the total number operator 
Tife. Since the IRHM commutes with the total number 
operator, _ff ^'^^ has the same eigenstates as IRHM! 

There are six open-loop hopping processes T^* depicted 
in Fig. [31 We analyze them sequentially below. 



E (1 - ^I^^) - 1 
E (1 - ^I^^) - 1 

k^i,l 

E(i-^I^^)-i 

k^i,l 

{N-X)-Y^b% 

3 

{N-1)-J2blb, 



b% 



E(i-^N 

{N-2)-Y^b% 
{N-1)-J2blb, 

3 

E(i-feK)-i 

k^i,l 

(^-2)-E^K 



b]b. 



b]b. 



b\b. 



b]b,. 



E Y.^M^'^^l^^ 
= E ii-bib,)Y.bib,b]b, 

= Y,(i-bib,)Y,b]b,b\h 



E(i-^I^^-) 
E%-i 



bjb. 



3 

{N-l)-Y,blbk) 



b\h.{m^) 



The hopping process Tj' in Fig. [3] (c) is expressed as 
= Ek^^.L, E,^u b%b]b4b, = Ti\ The fourth 
hopping process in Fig. [3] (d) is obtained as follows. 

= (i^b]b,)Y,blbkb\h 

j^ij.k k^i,l 



= Ti,' 



(B15) 



The hopping process Tg' in Fig. [3] (e) yields Tg* = 
T.]^^,Uk J2k^i,i blbjb\bkb]h = ri*. We analyze below the 
last hopping process in Fig. [21(f)- 

Tl'= E H^h^l^Abk 
= E ^I^'^ E ^l^M^^ 



E ^i^^ - 1 
E ^i^'^- - 1 

k^id 



E 

E%-i 



b]b. 



b\b,. (B16) 



(B13) 



We will now deal with closed-loop hopping processes 
such as those in Fig [H These lead to effective interac- 
tions. The process VI in Fig. [11(a), obtained from Fig. 
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[3] (a) by setting I = i, is given as follows. 



E(i-^I^^)-i 



E(i-^]^.) 



(^^-i)-E^I^* 



Next, the hopping process V2 corresponding to closed 
loop in Fig. |3] (b) is obtained from Fig. |3] (c) by taking 
I = i. 



E(i-^K) 



i 



(B18) 



Lastly, the hopping [depicted by the closed loop in 
Fig. m (c)] is obtained from Fig. |3] (e) by setting I = i. 

E T.^lbMbkb]b. 

= E(i-^^^^)E^M&. 



(B19) 



Finally, we consider Figs. [5] (a), (b), and (c) which deal 
with effective hopping terms Tp„ involving closed loops. 
The effective hopping term Tp\, corresponding to Fig. [5] 
(a), is obtained by setting fc = zin Fig. [3] (a): 



{N-2)-J2b]b, 



1^1 



{N-l)-Y^b]b, 



bjb. 



bjb,. 



(B20) 



To obtain the effective hopping term Tq^ corresponding 
to Fig. M (b), we take j = Hn Fig. |3] (a): 

7^C2 = E bWblbib]b, 

= ^ (1 - blbk)b\h 



k^l 



iN-l)-J2blb, 



b]b. 



b]h 



T'li 
ci- 



(B21) 

The effective hopping term T^g depicted in Fig. [5] (c) 
[upon setting k = i and j = / in Fig. [3] (a)] is given by 

Tc3 = b\b4bib\h = b\h. (B22) 
Thus we have shown that iJ*^^' contains effective hop- 
ping terms {Y.i.i>i T{Y,k'^k)b\bi + YI.C. ) and effective 
interaction terms (X^i ^(X]fc "-fc)"*)- Since T and V are 
functions of the total number operator, H'^'^^ and IRHM 
have the same eigenstates. These arguments can be ex- 
tended to even higher-order perturbation theory to show 
that the effective Hamiltonian (after taking all orders of 
perturbation into account) will give the same eigenstates 
as IRHM! 

We will now explain the expressions for the coeffi- 
cients tm Vn, and ten in Eq. (jBlip . obtained from 
third-order perturbation theory, using typical schematic 
diagrams shown in Fig. [S] [50|. We consider two dis- 
tinct time scales associated with hopping processes be- 
tween two sites: (i) ~ l/(Je~^ ) corresponding to ei- 
ther full distortion at a site to form a small polaronic 
potential well (of energy — g^wg) or full relaxation from 
the small polaronic distortion and (ii) ^ 1/J related 
to negligible distortion/relaxation at a site. The coef- 
ficient tn corresponds to the typical dominant distor- 
tion processes shown schematically in Fig. |6] (a) with 
the pertinent typical hopping processes being depicted 
in Fig. [3] (a). In Fig. M (a), after the HCB hops 
away from the initial site, the intermediate states have 
the same distortion as the initial state. Next, when the 
HCB hops to its final site there is a distortion at this 
final site with a concomitant relaxation at the initial 
site. Hence the contribution to the coefficient tn becomes 
J/{2g'^LOo) X J/(2g2cjo) x Je'S^ ~ J^e'^^ / {g^uj^f . As 
regards coefficient w„, it can be deduced based on the typ- 
ical dominant hopping-cum-distortion processes depicted 
in Fig. ini (b) which typifies the hopping processes in Fig. 
131(a). In Fig. HKb), when the particle hops to different 
sites and reaches finally the initial site, there is no change 
in distortion at any site. Hence Vn can be estimated to 
be J/(25^wo) X '//(S^^^o) x J ~ J^/(g^cjo)^. Lastly, we 
obtain the coefficient ten by considering the typical dom- 
inant diagram in Fig. [S] (c) corresponding to the typical 
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i 



v2/ 





j 

t 

j 





k 





k 





1 





1 



yOy • 



process in Fig. [S] (a). In Fig. IH] (c), where the first inter- 
mediate state depicts the particle hopping but leaving the 
distortion unchanged, we get a contribution J/(2(7^wo); 
for the next intermediate state, where the HCB returns 
to the initial site, the initial site has to undergo a slight 
relaxation (involving absorbing a phonon so as to yield a 
non-zero denominator in the perturbation theory) lead- 
ing to the contribution J/luq] and lastly, when the HCB 
hops to the final site, there is a distortion at the final site 
with a simultaneous relaxation at the initial site thereby 
producing a contribution Je~^ . Thus we calculate w„ to 
be J/{2g^uja) x J/wq x Jg-s" - J^e-a'' /{gujof [U. 





j 





k 



(a) 





k 





k 





j 





j 





k 



k 



j 



(b) 



^ ° ° 

< > I 



o • ° 

< J I 



^^ 

' J 1 



i J I 



(c) 

FIG. 6. Schematic diagrams (a), (b), and (c), corresponding 
to the hopping processes depicted in Fig. E] (a), Fig. |4](a), 
and Fig. [5] (a) respectively, yield coefficients t„, v„, and ten 
respectively. The intermediate states give the typical domi- 
nant contributions. Here empty circles correspond to empty 
sites, while filled circles indicate particle positions. Parabolic 
curve at a site depicts full distortion at that site with cor- 
responding energy —g^too {+g^Loo) if the hard-core-boson is 
present (absent) at that site. 
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